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ABSTRACT 

We formulate within a generalized distributional approach the treatment of the stability against radial pertur¬ 
bations for both neutral and charged stratified stars in Newtonian and Einstein’s gravity. We obtain from this 
approach the boundary conditions connecting two any phases within a star and underline its relevance for real¬ 
istic models of compact stars with phase transitions, owing to the modification of the star’s set of eigenmodes 
with respect to the continuous case. 

Subject headings: Thin-Shell Eormalism - Neutron Stars - Radial Perturbations — Phase Transitions 


1. INTRODUCTION 

There is theoretical evidence that compact stars, such 
as neutron stars, are made up of several matter phases 
(Shapiro & Teukolsky 1986). This is a consequence of the as¬ 
tonishingly high density excursion they could attain in their in¬ 
ner regions. Eor instance, a neutron star is in general thought 
as composed by at least to different regions: the crust and the 
core. Starting from very low densities of a few grams per cu¬ 
bic centimeter close to their surfaces, and up to densities of the 
order of the nuclear saturation value, pnuc ~ 2.7 x 10*"^ g cm“^, 
the crust of a neutron star is thought to be in a solid-like 
state. The core, with densities that might overcome by or¬ 
ders of magnitude pnuc, is instead thought to be in a liquid¬ 
like state. The details of the treatment of the thermodynamic 
transition (e.g. Maxwell or Gibbs phase construction), as well 
as the conditions on density and pressure at which such a 
transition occurs, are still a matter of debate. The applica¬ 
tion of the Gibbs construction with more than one conserved 
charge (e.g. baryon and electric) leads to the appearance of 
mixed phases, in between of the pure phases, with an equi¬ 
librium pressure that varies with the density, leading to a 
spatially extended phase-transition region of non-negligible 
thickness with respect to the star’s radius (Glendenning 1992; 
Glendenning & Pei 1995; Christiansen & Glendenning 1997; 
Glendenning & Schaffner-Bielich 1999; Christiansen et al. 
2000; Glendenning 2001). On the contrary, in the traditional 
Maxwell construction the phases are in “contact” each other. 
It is worth mentioning that in these treatments the pure phases 
are subjected to the condition of local charge neutrality, and 
so they do not account for the possible interior Coulomb fields. 
Indeed, the complete equilibrium of the multicomponent fluid 
in the cores of compact stars needs the presence of a Coulomb 
potential formed by electric charge separation due to gravito- 
polarization effects (Rotondo et al. 2011; Rueda et al. 201 1), 
favoring a sharp core-crust transition that ensures the global, 
but not the local, charge neutrality (see Belvedere et al. 2012, 
2014, and references therein). Besides the core-crust transi¬ 
tion, additional phase-transitions as the ones allowed by Quan- 
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turn Chromo-Dynamics (QCD) could occur within the core 
of the star itself (see, e.g., Glendenning 1996, and references 
therein). Erom all the above we can conclude that ultradense 
stars such as neutron stars shows necessarily a nontrivial strat¬ 
ification. Between any two phases, which can be very differ¬ 
ent, it is reasonable to investigate the situation where some 
quantities are discontinuous, e.g. the energy-density and the 
pressure. Such discontinuities can be harnessed by appropri¬ 
ate surface tensions. These surface quantities influence the 
stability of a system adding new boundary conditions to the 
problem which, as we shall show here, modify the set of eigen- 
frequencies and eigenmodes of a star. 

In this work we analyze the problem of perturbations in 
systems constituted of various phases that are split by sur¬ 
faces that host nontrivial degrees of freedom. This analysis 
is thought to be a generalization of the treatment for con¬ 
tinuous systems (see Herrera & Santos 1997, and references 
therein, for a comprehensive analysis of properties, types and 
stability of continuous anisotropic fluids also in presence of 
radiation and heat flux). By investigating the dynamics of 
perturbations, we are automatically probing the stability of 
systems. We shall trammel ourselves to the simplest possi¬ 
ble case: spherically symmetric extended bodies where radial 
perturbations take place. In order to model the problem, we 
shall assume these surfaces of discontinuity separating two ar¬ 
bitrary phases are very thin and a generalized distributional ap¬ 
proach (Poisson 2004; Raju 1982b) shall be adopted. We start 
our analysis in the Newtonian case in order to gain intuition 
of the relevant aspects of the problem, and finally we general¬ 
ize it to general relativity. Our purpose is solely to expound 
the problem and seek to solve it as generically as we can. Our 
analysis is far from being complete and has to be considered 
as a first step towards deeper investigations and scrutinies of 
specific cases shall be the object of posterior studies. 

We show in this work that phase transitions in the presence 
of surface degrees of freedom can be enclosed into additional 
boundary conditions to the problem. Besides, our formalism 
tells us that such boundary conditions are only self-consistent 
when the set of eigenfrequencies of the perturbation modes 
is related to the global system, not with individual phases. 
This is consistent with the well-known results from coupled 
springs, where there are only global frequencies. The pres¬ 
ence of further boundary conditions naturally modify the pos¬ 
sible set of eigenfrequencies, since we are inserting further 
restrictive aspects to the physical oscillation modes. There¬ 
fore, measurements on the pulsation modes in a star could tell 
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us very precisely about its internal structure, being a sort of 
fingerprint, which could help to understand better the nature 
of these systems. 

The article is organized as follows. In section 2 we formu¬ 
late the problem of radial perturbations of stars with interfaces 
in the classical Newtonian case. For the sake of exemplifica¬ 
tion, we compute the specific solution of the stability equa¬ 
tion in section 3 for the simple case of incompressible stars, 
scrutinizing which circumstances deliver us enough arbitrary 
constants to fix additional boundary conditions. In section 4 
we give the formulation in the case of systems with electric 
charge but still within Newtonian gravity. The formulation in 
the case of general relativity in the neutral case is presented in 
section 5, while in section 6 we extend the general relativistic 
treatment to systems endowed with electric charge. Finally, 
in section 7 we summarize and discuss the main results of 
this work. We use units such that c = G - 1 and the metric 
signature -2 throughout the article, unless otherwise stated. 


2. STABILITY OF CLASSICAL SYSTEMS WITH PHASE 
TRANSITIONS 

Assume a continuous classical astrophysical system with 
spherical symmetry. When its volume elements are perturbed 
radially, it is well-known (see, e.g., Shapiro & Teukolsky 
1986) that the evolution of perturbations of the form 

= (1) 


with w an arbitrary constant, is described by 


dr dr 


AdP , 
r dr 


( 2 ) 


where P{r) is the pressure of the background system under 
hydrostatic equilibrium and 


p dP 
P dp ’ 


( 3 ) 


with p{r) the mass density of the system. Formulated in this 
way, we have at hands an eigenvalue problem. For the case 
of continuous systems, the boundary conditions to be supple¬ 
mented to Eq. (2) are quite simple. They are directly related 
to the spherical symmetry of the system, as well as to the van¬ 
ishing of its pressure on its border even in the presence of per¬ 
turbations (other situations where a surface tension be present 
could also be envisaged and we shall attempt to elaborate on 
them in the sequel). In other words, we impose that 


^(0) = 0, and ^(Rj) = finite, (4) 


where Rs was defined as the radius of the star, such that 
P(Rs) - 0. For further details about these boundary con¬ 
ditions (Shapiro & Teukolsky 1986). Eq. (2) supplemented 
with Eq. (4) constitutes a Sturm-Liouville problem, where the 
aspects of its solutions are already known. Concerning the 
eigenfrequencies, they are all real and form a discrete hi¬ 
erarchical set. When one seeks for stable solutions to Eq. (2), 
one seeks for solutions with positive uP', specially for its fun¬ 
damental mode. As it can be seen from Eq. (1), negative val¬ 
ues of <jp indicate instabilities in the assumed background sys¬ 
tem, which leads to the conclusion that they do not linger on 
in time. They would either implode or explode. 

We now turn to the more involved problem of permitting 
the system to be stratified and harboring surface degrees of 
freedom on the interface of two given phases. Such degrees 


of freedom have themselves a dynamics, described generi- 
cally by the thin-shell formalism or Darmois-Israel formalism 
(Israel 1966, 1967; Lobo & Crawford 2005; Poisson 2004). 
We are here however particularly interested in another aspect 
of the problem, namely to understand the role such degrees of 
freedom play on the stability of the system when its parts (de¬ 
fined naturally by the hypersurface that hosts the aforesaid de¬ 
grees of freedom) are perturbed. Therefore, before anything, 
it is assumed that for being meaningful to talk about this sce¬ 
nario, one has that the hypersurfaces of discontinuity them¬ 
selves are stable (see Pereira et al. 2014, for further details). 
If this is not the case, any displacement of the hypersurface 
of discontinuity would trigger a cataclysmic set of events that 
would have as a result the disruption of the system. 

One expects that the stratified problem could be accounted 
for additional boundary conditions to the system. The reason 
for this is that the perturbations in the upper and lower regions 
with respect to a given surface of discontinuity would be de¬ 
scribed by the same physics [e.g. Eq. (2)], as well as the to¬ 
tality of the matches. The only missing points would be their 
connection (allowing combinations of solutions to be also so¬ 
lutions to the physical equations involved) and generalization 
by means of surface quantities. Eor example, the existence 
of a surface tension would account to an extra surface force 
term. The same ensues with the presence of a surface mass 
(enclosed by a surface mass density) and the associated pres¬ 
ence of a surface gravitational force. Therefore, in order to 
properly describe stratified systems, one must make use of 
distributions (Poisson 2004; Raju 1982b). 

We proceed now with the distributional generalization of 
the equations describing continuous fluids under gravitational 
fields. Assume that a surface harboring surface degrees of 
freedom of a system in equilibrium is at r = R. The first 
equation to be generalized in terms of distributions in this case 
is the equation of hydrostatic equilibrium. This should now 
read 

dP TP 

— +pg{r)-—d(,r-R)^^, (5) 

dr R 

where g{r) has been defined as the norm of the gravitational 
field, solution to (the distributional) Poisson’s equation V ■ g = 
—AnGp, g - -g(r)r, that can always be written as 

GM(r) r _ , _ 

g(r) = - 5 —, M(r) ± An I p(r)r^dr. (6) 

T Jo 

Besides, P in Eq. (5) stands for the surface tension (Peters 
2013) on the surface of discontinuity at r - R. The expres¬ 
sion on the aforesaid equation is the result of restoring surface 
forces (this is the reason they have an opposite direction to the 
pressure gradient) at a small surface area (the factor 2 comes 
from the principal curvatures in a surface element, which in 
the spherically symmetric case are equal). The gravitational 
force from the sheet of mass at r = R is naturally incorporated 
on the distributional definitions of p and g{r) as we shall show 
in the sequel. The mass density, on the other hand, must be 
expressed as 

p{r) - p~{r)9{R - r) + p^{r)0{r - R) + a6{r - R), (7) 

with 9(r - R) the Heaviside function, whose derivative is the 
Dirac delta function 6{r-R) in the sense of distributions. Prom 
the existence of a surface mass density, it can be checked that 
the gravitational field is discontinuous at a given surface of 
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discontinuity of the system (at R), i.e. 

[g(R)]l = 4nGcr, (8) 

where we have introduced the convention [A]l = - A , the 
jump of A(r) across r - R. Therefore, g{r) could be repre¬ 
sented distributionally as 

g = g^{r)e{r -R) + g~{r)e{R - r). (9) 

The above equation means that the associated distributional 
gravitational potential 0(r) (g - -V0, g = -gr) is always a 
continuous function, though not differentiable at a surface of 
discontinuity. 

Besides, we will assume that also the pressure can be dis¬ 
continuous at r - R, being hence written as 

F(r) = F-(r)ff(R - r) + F^(r)0(r - R). (10) 


The reason of why this is so will be clarified when we deal 
with our original problem in the scope of general relativity. 
The heuristic argument corroborating the validity of Eq. (10) 
is that it is meaningless to colligate a surface term to the radial 
pressure, since its associated force would necessarily be nor¬ 
mal to it and therefore would not lie on the surface. Only tan¬ 
gential pressures should have tied within their surface terms. 
From Eqs. (10), (7), (6) (8) and (9), we have 

and 

dF^ 

— +p^(r)g^(r)^0. ( 12 ) 

dr 

The arithmetic average present in Eq. (11) [see the def¬ 
inition of g{r) and the value of cr] is a general conse¬ 
quence of the product of delta functions with Heaviside 
ones in the generalized sense of distributions (Raju 1982b, a). 
Notice from Eq. (11) that its first term is the known 
Young-Laplace equation for spherical surfaces at equilib¬ 
rium (Rodrfguez-Valverde et al. 2003; Peters 2013, see, e.g.,), 
where just geometric aspects are taken into account for the sur¬ 
face tension. Its second term, though, is the gravitational sur¬ 
face tension, uniquely due to the non-zero surface mass. If the 
surface tension were null, the pressure jump could not be arbi¬ 
trary, but proportional to the surface mass density and must be 
a monotonically decreasing function of the radial coordinate 
[see Eqs. (8) and (11)]. From Eq. (11), one sees further that 
the force per unit area associated with the surface tension is 
exactly the one necessary to counterbalance both the forces 
coming from the pressure gradient at R and the surface gravi¬ 
tational force, as it should be. 

One sees that the above procedure generalizes our notion of 
hydrostatic equilibrium in each phase the stratified system has 
[see Eq. (12)] and automatically gives the surface tension at R 
that guarantees the hydrostatic equilibrium for arbitrary pres¬ 
sure jumps and surface masses. We will keep the same philos¬ 
ophy now concerning the generalization of Eq. (2). From our 
generalized hydrostatic equilibrium equation, we have that an 
important term for the deduction of the equation governing ra¬ 
dial perturbations would be the application of the Lagrangian 
operator A [AA = A(f, r + ^) - Ao(t, r), Aq and A being a 
physical quantity in the equilibrium and perturbed cases, re¬ 
spectively] on the surface force in Eq. (5), i.e. [see Eq. (20)] 


r 

j6ir-R) 


NP Pf 

±—5{r-R)--^5(r-R), (13) 


since IS.6{r - R) - 0 and AR = Now we assume that 
P - Pier). This means that we are endowing the fluid at 
the surface of discontinuity with adiabatic properties and the 
underlying microphysics is not contemplated in this proce¬ 
dure. For continuous media, the total mass in the interface of 
two phases is generally not a constant. This means that mass 
fluxes are allowed to take place. This generically would ren¬ 
der the mass of each phase not constant, an aspect not taken 
into account in Eq. (2). Nevertheless, if the displacements of 
the surface of discontinuity are small and oscillatory, we have 
that on average the masses on each phase are conserved (here 
it becomes clear why the surface of discontinuity should be 
stable). For adiabatic processes, we have 


AP = TfAcr, rf = -— , (14) 

ocr 

with 7]^ the square of the speed of the sound in the fluid at the 
surface of discontinuity. The missing term Act can be found 
via the thin-shell formalism when the classical limit is taken 
there. Generically, in the static and spherically symmetric 
case, cr can be written as (Lobo & Crawford 2005) 


4nGR 


[e-««>]!. 


(15) 


with the classical limit y6(r) a; GM{r)l(rc^) «c 1. It is easy to 
check that Eq. (15) reduces to Eq. (8) in the aforementioned 
limit. When perturbed, it can be shown that (3 ^ 13 + 6/3, with 
6/3 - -4TiGrpoe^^°^lc^ (Misner et al. 1973), po here mean¬ 
ing the mass density in the hydrostatic background solution. 
Hence, 

Act = 5cr-H O-qI =-[po|]^ + o-gf, (16) 

where we also considered ctq the background solution 
[Eq. (8)]. 

Another simpler way of obtaining Eq. (16) would be 
through the dynamics of 6g [6g= -(6g)r]. In the spherically 
symmetric case we have 6g - -47iGp^ (Shapiro & Teukolsky 
1986), and since 6g = g - go, with go the norm of the gravita¬ 
tional field without the perturbation f, from Eq. (8), we finally 
obtain 


AfT = + cr'o(R)i = <T'o(R)i - [po|]!. (17) 

47rO 

Besides, from Eqs. (6) and (8), for the case when the jump of 
^ is null at a surface of discontinuity (that shall be justified in 
the sequel), one shows that the above equation can be further 
simplified to 

Act = —cr|. (18) 


Now we show the general equation governing the propagation 
of radial perturbations. For F defined as in Eq. (10), p in terms 
of Eq. (7), g as given by Eq. (9), and finally 

|(r, t) = |+(r, t)e{r - R) + f9{R - r), (19) 


thus the equation governing the evolution of perturbations on 
a given volume element of the fluid is 


( dvr dF 


-^ + y(vv +v,) 


6ir-R)\ = 0. (20) 


Notice that we assumed that iv is a distribution like Physi¬ 
cally this must be taken because the phases are always “local- 
izable” and do not mix. In other words, this constraint reflects 
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the intuitive fact that the surface of discontinuity should be 
well-defined. The mathematical reason of why this is so shall 
be given below and it is related to the welTposedness of the 
problem. 

When developed, taking into account the hydrostatic equi¬ 
librium equation [see Eq. (5)], it can be simplified to 


and the boundary condition 

or (28) 

and therefore the only meaningful T are given by 

r(r) = r-(r)e{R -r) + r+(r)6»(r - R). (29) 


dr 


YP d 

yl Qy 




AdP~ 
r dr 



p (2,7V-3^) 


2 \ dfi dfi j 


6(r-R). 


( 21 ) 


First note that coefficients multiplied by 6^{r - R) or 6'(r - R) 
in Eq. (21) must be all null. Taking into account Eq. (19), it 
means that 


[|]! = 0. 


( 22 ) 


This automatically warrants ^ as a distribution without Dirac 
deltas, as we have advanced previously. In order to obtain 
Eq. (21), we used the results that for a distribution A{r) - 
A+{r)0(r -R)+ A-(r) 6(R- r). 


dA 

AA^6A + —i - [A]!f(5(r - R), 
dr 


\drj dr 
1 

^2 


(AA). 

d~C 


dj_ 
dr 


[A(R)]!b(r-R). 


(23) 


(24) 


Besides the above mathematical properties, we have also 
made use of 




dr 


2 \ dr 


dr 


(25) 


which is a direct consequence of assuming that the total mass 
in each phase is constant, even in the presence of perturba¬ 
tions. This is just guaranteed if the surface of discontinuity 
is stable, a prime hypothesis for having a well-posed stabil¬ 
ity problem. We also have assumed that P - P(p), which 
implies that = Y"^P"^j. In deriving Eq. (21), we 
further took into account Eq. (11). We finally stress that a 
simpler way to obtain Eq. (21) is to recall that AM = 0, 
which guarantees that Ag - —Ig^jr. The fact that AM - 0 
means that comoving observers with the fluid do not notice a 
mass change. The aforementioned result can also be directly 
showed by Eqs. ( 6 ), (24) and (25). 

It can be seen that only solutions of the type ^"^(r, t) - 
e™ 'f'^(r) for Eq. (21) are just meaningful if 

= oj. (26) 


This is the only way to eliminate the time dependence above 
inEq. ( 21 ), and also to guarantee that the jump of |be null for 
any surface of discontinuity at any time. Therefore, we arrive 
at the important conclusion that even a stratified system where 
oscillatory perturbations take place should be described by a 
sole set of frequencies. Each member of this set describes the 
eigenfrequency of whole system, instead of one or another 
phase. Nevertheless, we recall that at the surface of disconti¬ 
nuity the frequencies are in principle not defined. Bearing in 
mind the above conclusions, we have that, using Eq. (1), 

^(r) = r(r)0(R -r) + f^(r)0(r - R) (27) 


Gathering the above equations on Eq. (21), we obtain 


d 

dr 


YP d 
r^ dr 




AdP . 
r dr 




-m) 


^(3!P- 277V)-wV 


6{r-R). 


(30) 


Ones sees from Eq. (30) that in the case where P and cr are 
null, the classical expression, Eq. (2), is recovered. 

Summing up, substituting Eqs. (27) and (7) into Eq. (30), 
one sees that the only way to satisfy such an equation is by 
imposing that 


d 

dr 


1 d 
r^ dr' 


Y^P^-^{r^t) 


4 dP"'" , 7 + + 

r+ wV^^ = o, 

r dr 


(31) 


- 1 - 2{W - 277 V - 2P[P(P)]!)^(P) = 0, (32) 

and for completeness, condition (28). Eq. (31) is obtained 
here as a consequence of our distributional search for solu¬ 
tions to the radial Lagrangian displacements. This is exactly 
what one expects under physical arguments. Eqs. (32) and 
(28) are our desired boundary conditions to be further taken 
into account [besides Eq. (4)] at the interface of any two 
phases. 

For the case where [P(P)]1 = [r(f^)]t - cr - P - 0, we, 
have that also the derivative of ^ is continuous and therefore f 
is a differentiable function anywhere, as it should be since we 
are defining here a continuous system. Nevertheless, when¬ 
ever the aforementioned conditions do not take place, richer 
scenarios rise. Even in the case of a phase transition at con¬ 
stant pressure and negligible surface mass, the discontinuity 
of E and the existence of P generally render the derivative of 
f discontinuous. 


YP^ir^O 

dr 


3. A SPECIFIC EXAMPLE: UNIFORM DENSITY STARS 

We would like to stress that the boundary condition we have 
derived previously is actually very restrictive. The reason for 
this is that the physically acceptable cases [solutions to Eqs. 
(31) associated with a surface of discontinuity at r = R] are 
only the ones that deliver enough arbitrary constants of inte¬ 
gration for satisfying Eqs. (27) and (32). This should be taken 
into account together with the physical requirement of only 
admitting finite ^ everywhere and that are null at the origin 
[see Eq. (4)]. In the following we shall see a particular exam¬ 
ple where all of these aspects are evidenced. 

Let us now investigate a star made out of two phases, each 
with a uniform mass-density. Let us assume also, for the sake 
of simplicity and exemplification, that the associated E for 
each region is an arbitrary constant. As we will see, although 
this can be considered only as a first academic example, it al¬ 
ready evidences some aspects stratified systems should have. 
For this case it is straightforward to solve Poisson’s equa¬ 
tion and the equation of hydrostatic equilibrium (see, e.g., 
Shapiro &. Teukolsky 1986, for further details) and we have 
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for r < R 


P^{r) 


2nGp^_ 

~ 




3Po 

2nGp^ ’ 


(33) 


where is an arbitrary constant that corresponds to the pres¬ 
sure of system at the origin. For r > R, instead 


P+(r) 



(34) 


The constant mass-density in the inner and outer regions has 
been defined as p^ and p^, respectively. The pressure at the 
origin could always be chosen such that it matches the pres¬ 
sure at the base of the outer phase, and as it can be seen from 
Eq. (34), we have also introduced the condition of having a 
null pressure at the star’s surface. Substituting Eqs. (33) and 
(34) into Eq. (31), we are led to 


(1 -x\) 


dx\ 




(35) 


where we assumed that x+ = r/R,, x- = rjR., and 


. . 3^2 8 

“ 2nGp^T^ r± 


(36) 


We now solve Eq. (35) by the method of Erobenius. Eor a 
sake of simplicity, we shall drop the + notation. We therefore 
assume solutions of the form 


^ = (37) 

H=0 

where an and s are arbitrary constants to be fixed by primarily 
demanding that the first condition of Eq. (4) is satisfied, as 
well as f(x) is always finite. By substituting Eq. (37) into 
Eq. (35), it can be checked that the solutions to s are either 
s = 1 or s = -2. The associated recurrence relation obtained 
generally is 

{m +s){m +s+ 2>)-A 

^m+2 ~ f 

{m + s + 2)(m H- s -H 3) — 2 

with m - 0,2,4..., and a\ - aj, - - ... - 0. Let us analyze 

first the inner region. It is clear in this case that the associated 
flo for s - —2 must be null, as a consequence of one of our 
boundary conditions. Erom Eq. (38), one clearly sees that the 
power series given by Eq. (37) does not converge. Therefore, 
in order to satisfy the finiteness anywhere of we have to 
impose that the series be truncated somewhere, rendering it 
actually a polynomial. Hence 

A;,3=i =('«+!)('«+ 4). (39) 

Erom Eq. (36), one sees that just discrete frequencies [given 
by Eq. (39)] are possible to this region. Erom Eqs. (36) 
and (39), for having the frequency of the fundamental mode 
{m - 0) positive, one should have E^ > 4/3. Summing up, 
the physically relevant solution to this case just leaves out an 
arbitrary constant of integration, as required due to the scaling 
law present to ^ from Eq. (31). 

Let us now analyze the outer region. This is the most physi¬ 
cally interesting region since the problems at the star’s center 
are absent and therefore in principle one could even have two 
linearly independent solutions to Due to the finiteness of 


^ in this region, the outer counterpart of Eq. (39) must again 
take place. Nevertheless, for s - -2, one should also impose 

Kn.s=-2^im-2)im+l). (40) 

Erom Eq. (36), one sees from this case that its associated fun¬ 
damental mode {m - 0) is unstable. This means in principle 
that this solution to the outer region should be excluded, leav¬ 
ing out just the one from the case i = 1, where we should 
consider 

<.5=1 =('«+!)('«+ 4). (41) 

Notwithstanding, our previous analysis exhibits clear prob¬ 
lems; there are not enough arbitrary constants to fix Eqs. (28) 
and (32) and the eigenfrequencies in each region are differ¬ 
ent. However, we shall show that the condition of having a 
same eigenfrequency for the whole system, as required by our 
formalism, addresses all the problems. Obviously the stable 
eigenfrequencies of the star are only related to the solution 
i = 1. However, they could rise here either from aspects of 
the inner or the outer phases of the star. Let us see how this 
conclusion ensues. Assume initially that the only possible w’s 
are given by Eq. (39), associated with the modes . So, for 
having finite related to i = 1, one must impose that there 
exists a to the outer phase such that the numerator of the 
associated recurrence relation be null. It can be shown that 
this is just the case if 


-5+ ^9+AAUm:^J 


Therefore, Eq. (42) demands that 


(42) 


9 -H 4A;^i = {2p + \f, p>2, peN. (43) 

Eor the case s - -2 to it can be shown that the condition 
for the existence of a rsl^ted to fl,„+ = 0 is exactly 

given by Eq. (43). The mode itself is 


l+^9+4A:^,(m;^j) 


‘j=-2 


(44) 


Summarizing: if Eq. (43) is satisfied for any natural p > 2, 
there always exist modes, characterized by Eqs. (42) and (44), 
which guarantee the finiteness of as a linear combination of 
solutions for i = 1 and s = -2, associated with a given eigen¬ 
frequency a>in- ^ that just takes into account aspects of the in¬ 
ner phase of the system. In this case, one is able to come out 
with two arbitrary constants of integration, that would then 
guarantee that the additional boundary conditions raised by 
the stratification, Eqs. (28) and (32), be satisfied. It is im¬ 
mediate to see that a similar reasoning as above ensues if one 
chooses now w„,+ ^ as coming from aspects of the outer region, 
given now by Eqs. (36) and (41). Eor this case, we will find 
now a and a associated with w„,+j, as given by 

Eqs. (42) and (44), with the condition given by Eq. (43), re¬ 
placing A^^j(m/^j) by A/^j(ot^^j). Since and are given 
quantities, one sees that the only possible eigenfrequencies to 
system should satisfy 9 + 4A*^j = (2p -i- 1)“. This constraint 
is uniquely imposed due to the extra boundary conditions to 
the problem and is very restrictive. We have just shown a 
simple example where some of the aspects imprinted by strat¬ 
ification raise. Whenever there are two arbitrary solutions to 
^ in a given phase, it will be always possible to satisfy the 
constraints (28) and (32). 
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4. SYSTEMS WITH AN ELECTROMAGNETIC STRUCTURE 

Now we attempt to give a further step in our classical gen¬ 
eralization, by endowing the phases (as well as the surface 
of discontinuity) with an electromagnetic structure. Just for 
clarity, let us work with a system that exhibits just an electric 
field. The first point to be taken into account is the additional 
electric force present in the system. This would have the same 
structure as the gravitational force and therefore its generaliza¬ 
tion is straightforward. Now one should define also a distribu¬ 
tional solution to the charge density. The surface force asso¬ 
ciated with the surface tension should have the same form as 
previously, but now should also take into account the present 
electric aspects. The pressure in this case would also change 
due to the presence of the electric field and its jump over a 
surface of discontinuity could still be kept free. 

From the (distributional) Maxwell equations in the spheri¬ 
cally symmetric case one has that 

E{r) = Q{r) ± An { Pc{r)?^d?, (45) 

r- Jo 

where pdf) is the charge density at r. The associated “force 
density” is dFeildv - pcE{r)f. Therefore, the equation of 
hydrostatic equilibrium now reads 

dP Wn 

-j- + P(x)g{r) - pc{r)E{r) - -^6{r - /?) = 0. (46) 

dr R 

Therefore, like the gravitational field, the electric field also 
presents a jump at any surface of discontinuity (at r - R) 
endowed with surface charges. We write the charge-density 
as 


Pc{r) = p^ {r)e(R - r) 4 - {r)e{r - R) + (rc6{r - R), (47) 

while the distributional electric field is 

E{r)^ E^{r)6{R-r) + E^(r)6{r-R), [£(/?)]!= 4710-^. (48) 

By substituting now Eqs. (10), (45), (47) and (48) into 
Eq. (46), we have that the surface tension at equilibrium 
should read 

Po = ![/■(«)]; + (49) 

Notice that the existence of a surface mass would lead to 
[M^(R)]1 > 0, while [Q^{R)XL could in principle be any. The 
appearance of the last term in Eq. (49) is consistent with the 
expected and long ago known contribution of electric double¬ 
layers on the surface tension and surface energy of metals, 
as recalled by Erenkel (1917) in his seminal work. The exis¬ 
tence of such surface electric fields is well-known in material 
sciences and it has been determined experimentally from the 
photoelectric phenomenon by measuring the amount of work 
done by electrons to escape from the metal’s surface. There is 
a vast literature on the role of electric double-layers on surface 
phenomena in metals and contact surfaces, and we refer the 
reader for instance to Huang & Wyllie (1949); Israelachvili 
(2011), and references therein, for further details on this sub¬ 
ject. 

Since in the presence of an electric field the hydrostatic 
equilibrium equation and the surface tension changes, it can 
be checked that Eq. (30) keeps the same functional form. In 
drawing this conclusion, it was also assumed that the total 
charge of the system is a constant. This also means that 
AQ = 0. One also sees immediately that the main results 


concerning the stability of the stratified charged case are to¬ 
tally analogous to the neutral one, obtained by simply making 
the replacement P Pq. 


5. .STRATIFIED SYSTEMS IN GENERAL RELATIVITY 

Now we generalize the analysis of stratified systems to gen¬ 
eral relativity. Erom the classical analysis, we have learned 
that surface quantities must also be inserted into the gener¬ 
alized equation of hydrostatic equilibrium. Therefore, in a 
certain sense, we must find the proper generalization of the 
surface forces in general relativity. This will not be difficult 
bearing in mind the thin shell formalism, as we shall see in 
the sequel. Such a formalism states that in order to search for 
distributional solutions to general relativity, one has to con¬ 
sider an energy-momentum tensor at a surface of discontinu¬ 
ity, that we shall name S. It is precisely this surface content 
that leads to the jump of quantities that are related to physical 
observables, such as the extrinsic curvature. We now outline 
the formalism succinctly. Let us work just in the spherically 
symmetric case, where E is defined as = r - R{t) - 0, with 
T the the proper time of an observer on the aforesaid hypersur¬ 
face. Assume that the metrics in the regions above and below 
S (w.r.t. to the normal vector to it), described by the coordi¬ 
nate systems = (f±, r±, 6±, <p±), respectively, are given by 

dsl = - rld^l, (50) 

where 

= dOl + sin^ 6±d(p\. (51) 


Assume that the (three dimensional) hypersurface S be de¬ 
scribed by the (intrinsic) coordinates y° = (t, 0, ip) such that at 
the hypersurface f± = f±(T), 0± - 0 and ip± - ip, besides obvi¬ 
ously r± - R{t). In order to render the procedure consistent, 
one has to impose primarily that the intrinsic metric to X is 
unique. This fixes the coordinate transformations = ^(y^)- 
This is the generalization of the continuity of the gravitational 
potential across a surface harboring surface degrees of free¬ 
dom. Now, if the jump of the extrinsic curvature is non-null, it 
is automatically guaranteed the existence of a surface energy- 
momentum tensor (Poisson 2004) that in the spherically sym¬ 
metric case can always be cast as = diag(cr, -P, —P), with 
(see, e.g., Lobo & Crawford 2005) 


1 

AnR 




(52) 


cr 1 Ra'ie-^!^ + R-) + RR+B'RR^ 

P^ -+-- - - , (53) 

2 %nR yg-2)3 + ^2 J_ 

where generically A' = 5A/(9r and A = dAjdr. Einally, the 
discontinuity of the extrinsic curvature is the generalization 
of the discontinuity of the gravitational field across a surface 
with nontrivial degrees of freedom. The case of interest to 
be analyzed here is the static and stable (upon radial displace¬ 
ments of E) one R = R = 0, i.e. an equilibrium point. 

Let us see now how the generalization of the surface forces 
appear in this formalism. First of all, we know that 

T^y = T;^0{r -R) + T-^0(R - r) + e^ySat6{r - R), (54) 

where = dy‘‘ldx^ and Sab ^ hacSl, hat = e'^elgf.y. No¬ 
tice that is itself defined as a distribution. For hab, it does 
not matter the side of S one takes to evaluate it, since it must 
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be unique. Let us constrain ourselves first to the case of per¬ 
fect fluids (locally neutral) on each side of S. One sees from 
Eq. (54) and the coordinate transformations at S that 

± p = p+0(r -R)+ p-6(R -r) + o-6(r - R), (55) 

Tl ±-P = -P+6(r -R)- P eiR - r) (56) 

and 

tI = t] = -P, ^-P- r6{r - R). (57) 

From Eq. (56), we notice that there are not associated surface 
stresses. This is exactly what we advanced in the classical 
case with heuristic arguments and obtained here as a general 
consequence of distributional solutions to general relativity. 
Let us search formally for solutions to Einstein’s equations 
with the energy-momentum given by Eqs. (55)-(57) with the 
ansatz 

ds^ = - e^dp- - r^{de^ + sin^ Odip^). (58) 

The distributional nature of Eq. (58) will be evidenced by 
Eq. (54). As the solution, we have 

= 1-m(r) ^ 47r f p{r)T^df. (59) 

r Jo 

Notice from the above equation that p is given by Eq. (55) 
and therefore is a distribution generally discontinuous at 
R. For a we have, though 

a = ^[47rPr^ + m(r)], (60) 

ri 

where P is given by Eq. (56). From Eqs. (59) and (60), we see 
that a' is a distribution with no Dirac delta functions terms. 
Therefore, it implies that [a] 7 = 0. In other words, the 
function a is generally continuous though not differentiable 
at R. Nevertheless, from the conservation law of the energy- 
momentum tensor given by Eq. (54), we also have 

<(P± + EJ = -E1. (61) 

We notice that = 0, taking into account Eq. (54), would 
give us in principle terms dependent upon Heaviside func¬ 
tions, Dirac delta functions and their derivatives. The terms 
associated with the Heaviside functions are null due to the 
validity of Einstein’s equations on each side of E. The nul¬ 
lity of the remaining terms is associated with identities that 
the surface energy-momentum tensor has to satisfy (see, e.g., 
Mansouri (& Khorrami 1996). It is not difficult to show that 
such identities are automatically satisfied when one takes into 
account Eqs. (52) and (53) (see, e.g., Lobo (& Crawford 2005). 
This shows that the thin-shell formalism is consistent and the 
surface terms must indeed be taken as the aforesaid equations. 

In order to put Eq. (61) in the form that would allow us 
to consider Eqs. (55) and (56), we mandatorily should add 
surface terms. The correct way of doing it is 

dP 2'P t 

— = -(p + P)a' + —5(r - R) + \[P{R)]1 + 

cr cr \a'e^^Y\ 

- [<(R) + a'_{R)] + - - 4 ^^^~ }b(r - R). (62) 

Now we show that in the classical limit, Eq. (62) reduces ex¬ 
actly to Eq. (5) and thus it is its proper generalization. First of 
all, notice that in such a limit, a - 4>(r), 4>{r) the gravitational 


potential. Besides, it can be shown that in such a case 

where the above quantities are in cgs units. In Eq. (63) one 
recognizes the jump of the gravitational field g{r) - (p' - 
GM{r)lr^ at r - R, as exactly given by Eq. (8). Therefore, 


^{aUR) + aLm 


G[M^{R)]1 


(64) 

(65) 


Substituting Eq. (65) into Eq. (62), we see that the term inside 
the curly brackets of the latter equation is null [see Eqs. (53) 
and (11)]. Hence, the remaining term in front of the delta 
function is exactly 2'P jR as we already advanced and expected 
[see Eq. (5)]. 

Now we are in a position to talk about perturbations in the 
general relativistic scenario. When they take place, metric 
and fluid quantities change at a given spacetime point from 
their static counterparts. It is customary to assume that such 
departures are small, what allows us to work perturbatively. 
The primary task is to find such changes from the system of 
equations coming from relativistic hydrodynamics and gen¬ 
eral relativity. Nevertheless, these solutions are already very 
well-known (Misner et al. 1973). Our ultimate task is simply 
to generalize them to the distributional case. 

The equation governing the evolution of the fluid displace¬ 
ments on each side of S is the general relativistic Euler equa¬ 
tion, related to the orthogonal projection of (perfect flu¬ 
ids) onto i.e. 

(p -H P)u'',^u^‘ ±(p + P)a'' = (/’' - (66) 

where the labels “+" for each term in the above equation were 
omitted just not to overload the notation. 

In the hydrostatic case we have only that u'^ = e““o. When 
perturbations are present (Misner et al. 1973) 

=e-<(l-fc^), (67) 

where 6a is the change of the static solution ao in the pres¬ 
ence of perturbations at a given spacetime point. For the 
component, using the normalization condition u^ii^ - 1, one 
shows that (Misner et al. 1973) 

( 68 ) 

with f = d^'^ldP. Just for completeness, = 0. For 

the components of given by Eqs. (67) and (68), the left- 
hand side of Eq. (66) gives as the only nontrivial component 
a'j. - with 

- af - + 6a'^ + (69) 

and the associated equation of motion 

(p* + E^)(-flf) = -—. (70) 

dr± 

From Eq. (69), Eq. (70) can be cast as 

pp± 

(p^ + - (p^ + (71) 

Therefore, in terms of distributions Eq. (71) reads 
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dP 2'P 

(po + - IT - (P + - '^) + 

ar R 


[a'e-f^t 


cr 


“o*|+ + "o*!- + + Q''_| 


4;rR 
d{r-R), 


+ [PiR)]l 


(72) 


where p and P are given by Eqs. (55) and (56), respectively, 
and now 


|(r, t) = |+(r+, f+)6»(r+ - R) + ^(r-, r)6>(R - r'). (73) 


Notice that in Eq. (72) we are considering jumps and sym- 
metrizations of quantities defined in the presence of pertur¬ 
bations. As we stated previously, such perturbations change 
slightly the value of the physical quantities with respect to 
their hydrostatic values. The square brackets term of Eq. (72) 
is the proper generalization of the curly brackets term in 
Eq. (20). Naturally the reasoning for the eigenfrequencies 
of ^ in the general relativistic is the same as in the classical 
case. The same can be said about the continuity, though not 
differentiability, of f at any surface of discontinuity. 

Now, in order to have the proper generalization of Eq. (30), 
we should evaluate Eq. (72) at r+^ and then subtract it from its 
evaluation at r concerning the static solution. In order to do 
it properly, one should take into account the general results 
for the case of Lagrangian displacements coming from the 
standard procedure (see, e.g., Misner et al. 1973), but now in 
the sense of distributions, by recalling that AA = A{r + ^,t)- 
Ao(r, f), where Aq concerns the quantity A at equilibrium. 

It is not difficult to see that we have the following re¬ 
sults in the general relativistic distributional case (see, e.g., 
Misner et al. 1973, for the treatment of a continuous system): 




g2/3o 

Aa' =47rre2^“ [6P + IdpPo] + - Sp + a'J, 

r 


AP^-yPo 



+ 6P 


Ap = -(po -I- Po) 


1 d 


\ 2 I dr) 


(74) 

(75) 

(76) 


+A(t+^^ + [Po]!|)b(r - P) - 
. po + PQ (dP\ 

^ Po Up 


Acr = - 


2cro| 

R 


-I 



-H 


[coshjSo]!) 

47rP2 j ■ 


(77) 

(78) 

(79) 


We stress that Eq. (76) assumes adiabatic processes, in which 
one considers P = P(p) and we have made use of [^]t - 0 for 
the above equations. 

We shall seek for solutions to the perturbations as | = 
e“^‘^(r) with to the same for all the phases the system may 
have. We just need to worry about the Dirac delta function 
term, since it gives us the desired boundary condition valid 
for the separation of each two phases. The terms in front of 
the Heaviside functions by default shall be the ones found in 
continuous media. It is not hard to see that the surface terms 
at the end should satisfy the condition 


Act , A[a'e^]^ 

-(277“ -Hi) - — - - = 0, 

R^ ' ^ 47tR 


(80) 


where we have used Eq. (24). When the last term on the left- 
hand side of above equation is expanded, by using Eqs. (60) 
and (74), Eq. (80) can be further simplified to 


In^Ao- , 
A 

[coshyS]); 


R 

AnR^ 







[APe^-]!. 


( 81 ) 

One sees from Eq. (81) that Eq. (32) is recovered in the clas¬ 
sical limit by recalling that p - M{r)lr, which implies that 
the last term on the left-hand side of the above equation is 
4cr{g^ + g^)l{2R). Besides, in this limit we take P —> 0 and 
^ 1 for the remaining terms. 

The case where electromagnetic interactions are also 
present is also of interest since its associated energy- 
momentum tensor is anisotropic. This naturally influences the 
equation of hydrostatic equilibrium, since it now becomes 


a'(P-Hp) = -P'- 


2 

-iP, - P), 

r 


(82) 


where P, P, and p are the resultant radial pressure, tangen¬ 
tial pressure and energy density of the fluid, respectively. Eor 
the electromagnetic fields, clearly (P, - P) is solely related 
to them. Due to the aforementioned aspects, the latter should 
also influence the dynamics of the radial perturbations, as we 
shall show in the next section. 


6. ELECTROMAGNETIC INTERACTIONS IN STRATIFIED SYSTEMS 
WITHIN GR 


We consider now the inclusion of electromagnetic interac¬ 
tions within the scope of stratified systems in general relativ¬ 
ity. An important comment at this level is in order. Since we 
are dealing with electromagnetic fields in stars, it would be 
more reasonable to assume the Maxwell equations in material 
media. Nevertheless, since the knowledge of the structures 
constituting the stars is not yet precise, it is difficult to as¬ 
sess their realistic dielectric properties. Since working with 
Maxwell equations in the absence of material media gives 
us upper limits to the fields under normal circumstances, this 
seems to be a good first tool to evaluate the relevance and ef¬ 
fects of electromagnetism in stars. We shall adopt for the time 
being to follow this approach. The energy-momentum tensor 
of each layer of the system we are now interested should also 
have the electromagnetic one i.e. 

4777 ^ 7 ) = -F,aFypg‘^‘' + g,y , (83) 


where we defined P^y = dfjAy - dyA^^. Solving Einstein- 
Maxwell equations on each layer of a stratified system leads 
us to the following equilibrium condition (Bekenstein 1971) 


dP 


Q{r)Q'{r) 

47rr^ 


- a'gip + P), 


(84) 


where 

Q(r)± f Anr^p.Adr, £(r) = e“e+- 8 e.^^ ( 35 ) 

Jo 

^ We restrict our analyses to the Maxwell Lagrangian, - -F/4. 
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Anr^P + mQ(r) — 


Q\r) 


and 


r ? 1 r 

iQir) - I 4nr^pdr +^ + :^ —dr. 
Jo 2r 2 Jo r- 


( 86 ) 

(87) 

( 88 ) 


Eqs. (85), (86), (87) and (88) are the charge, the radial and 
time components of the metric and the energy of the system up 
to a radial coordinate r, respectively. Besides, in Eq. (85), we 
also exhibited the electric field E(r) in the context of general 
relativity, obtained by means of the definition Ftr - E{r) - 
-drAo. We stress that pc is the physical charge density of the 
system, defined in terms of the four-current by 
where is the four-velocity of the fluid with respect to the 
coordinate system (f, r, 6, ip) (see Landau & Lifshitz 1975, for 
further details). 

Erom Eq. (84) one sees that in the scope of general relativ¬ 
ity the effect of the charge is not merely to counterbalance the 
gravitational pull. Eor certain cases, it could even contribute 
to it. The reason for that is due to the contribution of the elec¬ 
tromagnetic energy to the final mass of the system, as clearly 
given by Eq. (88). Notice from Eq. (88) that we have assumed 
that the mass at the origin is null, in order to avoid singular¬ 
ities there. More generically, one could assume point or sur¬ 
face mass contributions in Eq. (88) by conveniently adding 
Dirac delta functions in p. Einally, we stress that Eq. (85) 
can indeed be seen as the generalization of the charge in gen¬ 
eral relativity, since it takes into account the nontrivial con¬ 
tribution coming from the spacetime warp due to its energy- 
momentum content. 


Notice that the classical limit to Eq. (84) can be shown to 
coincide with Eq. (46), by recalling that EdaJr) — Qcias{r)lr^ 
and, from Eq. (85), - Anr-pc- Besides, we recall 

that when brought to cgs units, the term Ir (here in geo¬ 
metric units) becomes Hc^r), which is null in the classic 
non-relativistic limit, as well as any pressure term on the right- 
hand side of the aforementioned equation. 

Now, consider the analysis of a charged system consti¬ 
tuted of two parts, connected by a surface of discontinuity (at 
- R) which hosts surface degrees of freedom, such as an 
energy density, charge density and a surface tension. Its gen¬ 
eralization to an arbitrary number of layers is immediate since 
each surface of discontinuity is only split by two phases. The 
proper description of the charge density in this case would be 
given by the generalization of Eq. (47). Therefore, one would 
have at equilibrium that 

Q{r) = Q-{r-)e{R - r) + e+(r+)0(r+ - R). (89) 

and for Q'{r), a Dirac delta shall rise, due to pc- 

Let us define the distribution 

pc = pce^° - ptS(r - R)+ PcO(R - r) + d-c6(r - R), (90) 

where p^- p*(r=^). Erom the above definition, we have that 
Q' = Anr^pc- It implies that the total charge is the same as 
the one associated with pc in an Euclidean space. Therefore, 
all classical results apropos of the charge densities and total 
charges that we deduced in the previous sections ensue here 

for Pc- 

We seek now for the distributional generalization of 
Eq. (84). This can be easily done by following the same rea¬ 
soning from the previous section, which finally leads us to 


dP 

dr 


^^-(p + P)a'Q+^6ir-R)+l [P{R)]: 


R 




CTq CTq 

— [a'{R) + a' (R)] + — - - - 


^_[Qm + Q-(R)]'j6{r-R), 
(91) 


where we are assuming surface quantities with the subindex 
“2” are related to the charged versions of Eqs. (52) and (53) 
[see also Eqs. (86) and (87)]. It is easy to show that in the 
classical limit Eq. (49) naturally raises, implying that in such 
a limit the curly brackets in Eq. (91) is null. 

We consider now the case where radial perturbations take 


place in our charged system. This case is more involved 
than the neutral case since the charged particles also feel 
an electric force. The equation describing the evolution of 
the displacements can be shown to be generalized to (see 
Anninos & Rothman 2002, for the dynamics of the radial per¬ 
turbations in a given phase) 


dP 




Q{r)pc 

H- —o(r - R) ■ 


R 


R 


AnR 


■ “2+ + “e- 


^ [er(-R) + e-(-R)] + [P(R)]- 6{r - R). 


(92) 


Eor the change in Q{r), it can be shown (see Bekenstein 197 1) 
that in the comoving frame there are no currents. This means 
that the Lagrangian displacements of Q are null, AQ - 0. 
Eq. (92) takes into account the values of the physical quan¬ 
tities in the presence of perturbations at r. In order to obtain 
the generalization of Eq. (30), we should evaluate Eq. (92) at 


r + ^ and subtract it from Eq. (91). This is due to the definition 
of the Lagrangian displacement of a given physical quantity, 
intrinsically related to the notion of comoving observers with 
the fluid, who naturally could describe its thermodynamics. 

In order to simplify Eq. (92), we have that in the gener¬ 
alized charged case (see Anninos & Rothman 2002, for the 
treatment in a phase of a charged system) 



















10 


Pereira and Rueda 




AtTg = 47rre^^2" 


(93) 


6P + ISPq 


ri) — 


AP =-yPo 


-Poa 


(r^e^eoiy 


Ap = -(po + Pq) 
po + Po (dP 


1 d 


STrr'* 
+ 6Pq 


5 47re2^eop^2o j 27ro-c| 

- ^Pq + f “ 


R 


^ ^ 

dr r^ 


I^QOi 

idr 

\ 2 ' 

[ dr 


r ^ 


Acto = - 


Po \5p 

2cr go 


R 


[g/'eopo]) 


[coshjSgoll 


47rR2 




(eoe'^»+ eo^’'^»>(^-p), (94) 

(95) 

• [PoT-l - §1 [ 2 + + 2 -] ) 5 (^ - R), 

(96) 

(97) 

(98) 


R 


The additional “0” subindex in a physical quantity means that 
its value at equilibrium was taken. We just stress that Eq. (97) 
is the general relativistic definition of the adiabatic index and 
assumes the existence of an equation of state linking the pres¬ 
sure and the density of the system, P - P(p). In this sense, it 
generalizes E as defined by Eg. (3). 

By seeking for solutions ^ = e'"'f(r), one can see that 
Eq. (92) just gives meaningful boundary conditions when a 
frequency a> is the same for all the phases present in the sys¬ 
tem. We emphasize this is an universal property of the ap¬ 
proach developed here, due to the surface degrees of freedom 
and the well-posedness of the problem of radial perturbations 
in stratified systems. The associated boundary condition ris¬ 
ing from this analysis leads us to the conclusion that genet¬ 
ically ^(r) is not differentiable at a surface of discontinuity, 
though continuous. It can be shown that the associated bound¬ 
ary condition to be taken into account here is functionally the 
same as Eq. (80) [or Eq. (81)], where now the metric and sur¬ 
face quantities should be related to the charged case. 

7. CONCLUSIONS 

In this article we have developed a formalism for assess¬ 
ing the stability of an stratified star against radial perturba¬ 
tions. We have derived the relevant equations defining this 
boundary-value problem for both neutral and charged stars 
and in Newtonian and Einstein’s gravity. It makes use of the 
generalized theory of distributions, since we assumed that the 
surfaces of discontinuity are thin, that they host surface de¬ 
grees of freedom, and that the phases separated by them do 
not mix. We showed that though the phases may be very differ¬ 
ent amongst themselves, when perturbations take place, they 
lead to the notion of a set of eigenfrequencies describing the 
whole system, instead of an independent set for each phase. 
Besides, our formalism gave us as a consequence the proper 
additional boundary conditions to be taken into account when 
working with stratified systems. Such boundary conditions 
encompass surface degrees of freedom in the surfaces of dis¬ 
continuity and generically modify the set of eigenfrequencies 
with respect to their continuous counterpart. This should be 
a generic fingerprint of stratified systems with nontrivial sur¬ 
face degrees of freedom. Our analyses are relevant for the 
assessment of the stability of realistic star models, since they 
ensue the precise notion of boundary conditions. It was not 


our objective to systematically apply our formalism here, but 
simply to derive and expound it. It is clear that for precise and 
realistic numerical stability calculations, it would be ideal to 
have the microphysical knowledge of the properties of the in¬ 
terfacial surfaces. However, this a difficult problem which has 
been elusive even in the most advanced field of material sci¬ 
ence, where laboratory data of the material surface properties 
are accessible, but still with the lack of a complete physical 
theory for their explanation (Israelachvili 2011). Thus, the 
measurement of the star’s eigenmodes becomes of major rel¬ 
evance, since it could give information not only on the star’s 
bulk structure but also on the possible existence of interior 
interfaces and their associate microphysical and electromag¬ 
netic phenomena. 

The radial instabilities shown in our analyses should be in¬ 
terpreted analogously as for continuous stars since, even in the 
stratified case, a global set of eigenfrequencies raises. Thus, 
stratified stars would either implode or explode when they are 
radially unstable. Though we were concerned only with radial 
perturbations, it is of interest to investigate the additional os¬ 
cillation modes owing to nonradial perturbations. The only 
point to be added, with respect to continuous stars, is the 
proper redefinition of the surfaces of discontinuity when such 
perturbations take place. This is clearly a richer scenario that 
insert additional degrees of freedom into the system, leading 
to the appearance of additional modes such as the gravita¬ 
tional ‘g-modes’ (see, e.g., Reisenegger (fe Goldreich 1992). 
Such an analysis, however, is as a second step that goes be¬ 
yond the goal of the present work, and that we are planning to 
investigate elsewhere. 
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